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1. INTRODUCTION 


The screwworm is the larval stage of the fly Coohliomyig hnrninivorax , Coquerel . 
The larvae feed only on living tissue forming a parasitic relation with a 
wide range of wild and domestic animals. Tick bites, wounds left by castra- 
tion anu dehorning, and the navals of newborn animals are common points of 
Infestation. The larvae move in place with a corkscrew motion, hence the 
name screwworm. Left untreated, crippling or death of the host animal is 
Inevitable. 

Screwworm flies at one time were found from Florida to California and as far 
north as Nebraska. In the late 1800' s, many lifestock producers were occa- 
sionally forced to restock due to heavy losses caused by the screwworm. 

Before effective controls, domestic livestock losses were estimated at 
$20 million _...iually In the Southeastern United States and $50 to 
$100 million in the Southwest. 

In the early 1 950 ' s , researchers with the U.S. Department of Agriculture 
(USDA) showed that careful irradiation of the pupal stage of the fly produces 
sexual sterility (ref. 1). The sterile flies continue to mature normally and 
demonstrate the same types of activity as wild flies. It was noted by 
researchers that the typical female screwworm mates only once, while the male 
continues to mate throughout its lifetime. When large numbers of sterile 
screwworms are introduced into the environment, a sizable portion of the wild 
female flies mate with sterile males resulting in infertile eggs. Thus, the 
fly's own reproductive cycle is used to limit its populations. 

Efforts to completely eradicate screwworm populations through animal husbandry 
practices seem to have failed because wild animal populations maintain a 
screwworm population which can reinfest domestic animals. In 1955, it was 
estimated that 50 to 90 percent of the screwworm flies in Texas were produced 
by wild animals. More recent estimates place this fiqure much lower. Even 
so, if only 20 percent of an existing screwworm population reproduce, the 
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next generation will be the same size as the parent generation. Thus, a 
program such as the USDA sterile male drops, which affect the entire repro- 
ductive population. Is the only approach which can hope to eradicate the 
screwworm. The sterile-male program has clearly resulted In a sharp reduction 
In the number of screwworm Infestations In both the United States ard Mexico 
(ref. 2), yet a number of puzzling outbreaks have occurred. 

An excellent survey of the research leading to the USDA program for screwworm 
eradication Is given In reference 3. More detailed Information may be found 
In reference: 4 and 5. 

In 1973, the Health Applications Office at The National Aeronautics and Space 
Administration, Lyndon B. Johnson Space Center (NASA/JSC) began a program to 
evaluate the use of remotely sensed data as an additional tool In existing 
and projected eradication efforts. The need tc relate remotely sensed data 
to screwworm Infestations resulted in a two-part developmental study. First, 
the ability to use remotely sensed data to estimate weather conditions was 
evaluated. Second, the effect of weather on screwworm populations was modeled. 
Barnes and Forsberg (ref. 6) have provided an excellent overview of the pro- 
gram's approach. This report deals with the salient points of the weather- 
population interaction. 
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2. ROLE OF WEATHER IN SCREWWORM GROWTH AND SURVIVAL 


That weather conditions, particularly soil moisture and air temperature. 
Influence the growth and survival of screwworm populations has been recog- 
nized for a number of years (ref. 7). Meteorological extremes have served to 
limit the distribution of S' rewworm infestations. In fact, the unusually 
cold winter of 1957-58 has been credited with a large contributory role In 
the successful effort to eradicate screwworms from the Southeast (ref. 8). 
Currently, the screwworm is restricted to tropical and subtropical areas of 
North and South America during the winter months, with large populations over- 
wintering each year in Mexico. Each summer, as temperatures Increase, screw- 
worm flies advance from Mexico into the United States and remain until cool 
temperatures prevail. Figure 2-1 shows the total number of screwworm infesta- 
tions reported in Texas, by month, since the beginning of the eradication 
program in 1962. Of particular interest are the outbreaks occurring in 1968 
and 1972. These outbreaks have, at least in part, been attributed to favor- 
able weather conditions. 

As a point of departure, an extensive literature review of the effect of 
selected environmental conditions on the growth and development of the screw- 
worm was undertaken (ref. 9). Based on that study, it was recognized that a 
wide range of meteorological parameters can affect fly activity including 
air and soil temperatures, precipitation, sky cover, humidity, windspeed, and 
soil moisture. In addition to meteorological conditions, there are a number 
of other important factors which influence fly populations: availability of 

hosts, predators, incidence of wounds, topography, condition and type of 
vegetation, and the effects of the eradication program itself. 

All of these concomitant factors can and do play a role in the observed 
fluctuations of screwworm populations. However, until recently, there has 
been no real way to quantify the effect of a given set of weather conditions 
on screwworm populations. Historically, attempts to relate weather conditions 
to the number of screwworm infestations, while producing mixed results, have 
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clearly Identified air temperature and soil moisture conditions as Important 
parameters (ref. 10). The problem then Is to quantify the effects of a 
given set of meteorological conditions such that standard statistical tech- 
niques can be used to analyze weather-screwworm relationships. 

2 • 1 SOIL M OISTURE AND AIR TEM PERA TURE 

As might be expected, the nature of the available data shaped the basic study 
of weather-screwworm relationships. Data on several of the desirable param- 
eters were either completely unavailable or available only at prohibitive 
cost. 

Based upon the literature survey and upon a number of preliminary statistical 
studies. It was decided to concentrate on the effects of soil moisture and 
air temperature. 

Soi 1 -moisture data based upon dire- t measurement are extremely limited and 
generally unavailable. The use of a moisture-budget approach, utilizing 
more commonly measured meteorological variables, has been suggested by a 
number of authors. The approach using conventional measurements of air 
temperature and precipitation (ref. 11) was extended to provide a measure of 
long-term moisture anomalies (ref. 12). 

In 1968, a related measure of short-term moisture anomalies was also developed 
(ref. 12). This parameter, called the Crop Moisture Index (CMI), has come 
into general usage and is published weekly during the summer in the Weekly 
Weather and Crop Bulletin of the U.S. Department of Commerce. The CMI varies 
around zero. When positive, the index indicates a surplus of moisture and 
conversely a deficit when negative. 

The CMI is calculated on a weekly basis for each climatological division in 
the United States. These divisions are multicounty areas considered to be 
fairly homogeneous both cl imatological ly and hydrological ly. The 10 clima- 
tological divisions for Texas are shown in figure 2-2. 
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Figure 2-2.- Texas cl imatogical divisions. 


The meteorological data base for this study consisted 01 weekly values of 
mean air temperature, total precipitation, and CM I for each climatological 
division In the Southwestern United Statist. Tnese values represent the 
average of data from all meteorological stations in the division. 

The basic Infestation data were taken from weekly Agricultural Research Science 
summar<es by county for the years 1967-73. The data were further summarized 
Into weekly to.al Infestations for each climatological division. 

The study of environmental influence on the screwworm Is complicated by 
several factors. First, data on the number of cases are based primarily on 
rancher reports of infected livestock. These reports may be delayed for 
various reasons, and some cases may never be reported. When and how often 
cases "occur" Is partly a function of how frequently a rancher checks his 
herds, a factor controlled partly by previously reported casts in his area. 
Thus, a type of "feedback" exists, where a greater percentage of actual cases 
is likely to be reported during a serious outbreak than during a time when 
screwworm activity is low. 

An additional complication is the eradication program itself, which holds 
uown the amount of screwworm activity that might otherwise occur under a 
given set of environmental conditions. Because of the nature of the eradi- 
cation program, it is not possible to completely separate the variation due 
to program from that due to environmental parameters. 

2.2 METHODOLOGY DEVELOPED TO DETERMINE WFATHER INFLUENCE ON SCREWWORM GROWTH 
AND SURVIVAL 

Based upon data obtained from the literature survey and the screwworm's life 
cycle as shown in figure 2-3, a sample model was developed to examine the 
influence of weather on screwworm growth and survival. This model was used to 
calculate the average weather and resulting growth per generation for several 
years of historical data in the Southwestern United States. 
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Figure 2-3.- Idealized diagram of the screwworm life cycle. 


The growth per generation was defined as the ratio of the number of cases 
during a given week to the number of cases reported during the week In which 
the parent generation occurred. The length of time between the parent and 
progeny generation was determined using the average times given for each phase 
In figure 2-3. The length of the pupal phase, which Is temperature dependent, 
was determined by a degree-day calculation. 

Laboratory data on the development time of screwworm pupae at various temper- 
atures are ..hown In figure 2-4. By converting these developmental times to 
their reciprocals and extending the resulting line to zero, a developmental 
threshold value of 10.8° C (51.44® F) was estimated along with a thermal 
constant of 139.1 degree-days. The length, n, of the pupae phase Is the 
first n which satisfies the following equation 

£ (T 1 - 10.8) s 139.1 
L-l 

where T^ Is the daily mean air temperature In °C. 

Average conditions of temperature, total precipitation, and CM I value were 
computed for the adult phase prior to first ovlposlti' n and the pupal phase. 

The weather conditions during each phase we^e then plotted against the result- 
ing growth per generation. 

The plots of weather conditions versus screwworm population growth resulted 
In the definition of a set of screwworm potential functions. This was 
accomplished by observing the maximum growth per generation under a wide 
range of temperature and soil -moisture conditions. For conmonly repeated 
conditions, a wide range of growths was observed. Those growths, less than 
the potential, can be attributed to a variety of limiting factors including 
a lack of wounds, natural predators, and effects of the eradication program 
itself. Some points were observed which represented unexpectedly large 
growths. These points were attributed to sampling error In some cases and 
to the migration of flies Into areas to take advantage of favorable conditions. 
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Figure 2-4.— Laboratory data on the development time of screwworm pupae at various temperatures. 



Careful analysis of all available data resulted In the weather potential 
functions given In figures 2-5 and 2-6. These figures are Intended to repre- 
sent the potential growth of a single generation of screwworm, If that meteoro- 
logical parameter Is the only factor acting to limit the screwworm population. 

The screwwtrm potential, given as a function of air temperature, has a gradual 
increase from zero near 10° C (50° F) to a maximum of approximately 27° C 
(82° F), followed by a sharp decrease around 30° C (85° F). 

The screwworm potential as a function of soil moisture Increases to a maximum 
near 0.8 anl th> .i falls with Increasing moisture. This function represents 
an empirical verification of laboratory reports showing a lower survival for 
screwworm pupae under extremes of soil -moisture conditions. 

In a similar manner, the screwworm potential as a function of generation 
length was derived. Shown In figure 2-7, the potential Is maximum when the 
generation length is near 3-1/2 weeks. 

To evaluate the potential utility of these growth functions four models were 
developed. Two model- were designed to predict growth per generation, 
whereas the remaining models predicted actual infestations. Each model was 
built using the same set of data used to develop the growth functions. 

The first model, known as Model 1, is a multiplicative model for growth 
and the formula follows: 


d l d 2 a 3 a 4 
G = ax, x 9 x- J x. 
o ol234 

where 

G o - growth per generation 

Xj - growth potential of mean air temperature during the adult phases 
*2 — growth potential of mean air temperature during the pupae phases 
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MEAN SOIL MOISTURE A. MEASURED BY 
THE PALMER CROP MOISTURE INDEX 


RELATIONSHIP BETWEEN MEAN AIR TEMPERATURE 
AND GROWTH PER GENERATION 
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Figurp 2-6.- Relationship between mean air temperature and growth per generation. 



m 



IVliNBlOd H1MOH0 


2-12 


GENERATION LENGTH (DAYS) 

Figure 2-7.- Relationship between generation length and growth per generation. 



x- — growth potential of CMI 

«9 

x.- growth potential of mean pupation time 

The second model Includes the Infestation, as measured by reported cases, of 
the parent generation as an additional Independent variable. Model 2 may be 
written: 


b x b l b 2 b 3 b 4 b 5 
0 1 x 2 x 3 x 4 x 5 


where 


I -current infestation 
o 

x,. — infestation during parent generation 


The third model adds the actual growth per generation observed 1 week in the 
past. Model 3 may be written: 


c X^xVVVC 6 
C 0 X 1 X 1 x 2 x 4 x 5 x 6 


where 


x & - observed growth per generation 1 week in the past 


The final model evaluated utility of Model 3 to predict infestation. 
Model 4 may be written: 


Each of these models is either linear or implicitly linear. Thus, the models 
may be fitted through multiple linear regression. Tables 2-1 through 2-4 
present a summary of statistics for each model build. 

Model 1 predicts the potential growth generation as a function of observed 
weather conditions. The adjusted multiple correlation coefficient (R) was 
0.515 when the model was fitted to the historical data set. The growth 
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TABLE 2-1.- STATISTICS FOR MODEL 1 
(a) Regression coefficient 


Variable 

Coefficients* 

Standard error 

t-value 

X 1 

0.32317 

0.19940 

1.62 

x 2 

-.17854 

.27203 

-.66 

x 3 

x 4 

Jin (constant) 

.37625 

.12131 

3.10 + 

.57293 

-71.06738 

.19430 

32.67852 

2.94" 


(b) Analysis of variance 


Source 

Degrees of freedom 

Sum of squares 

Mean of squares 

F-value 

Regression 

4 

233613.09 

58403.27 

16.34 

Res i . jal 

166 

593292.92 

3574.06 


Total 

170 

826906.01 




Significant at 0.95 level. 

★ 

Adjusted multiple correlation coefficient: 0.515. 
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TABLE 2-2.- STATISTICS FOR MODEL 2 
(a) Regression coefficient 


Variable 

Coefficients* 

Standard error 

t-value 

X 1 

x 2 

x 3 

0.15213 

0.21645 

0.70 

-.21072 

.53777 

.29138 

.13836 

-.72 

3.89' 

+ 

x 4 

.79839 

.21263 

3.75 

x 5 

.70187 

.06839 

10.26 + 

en (constant) 

.43349 

.29920 



(b) Analysis of variance 


Source 

Degrees of freedom 

Sum of squares 

Mean of squares 

F- value 

Regression 

5 

68.92 

13.78 

34.70 

Residual 

163 

64.75 

.40 


Total 

168 

133.67 




Significant at 0.95 level. 

Adjusted multiple correlation coefficient: 0.708. 
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TABLE 2-3.- STATISTICS FOR MODEL 3 
(a) Regression coefficient 


Variable 

Coefficients* 

Standard error 

t-value 

X 1 

x 2 

x 3 

x 4 

x 5 

x 6 

in (constant) 

0.01236 

-.16318 

.23395 

-.48307 

-.20947 

-.59515 

.432.14 

0.09462 

.21757 

.12736 

.19035 

.06269 

.07180 

.28467 

0.13 
-.75 
1.84 + 
2.54 + 
-3. 34 + 
8.29 


(b) Analysis of variance 


Source 

Degrees of freedom 

Sum of squares 

Mean of squares 

— 

F-value 

Regression 

6 

34.86 

5.81 

27.24 

Residual 

117 

24.95 

.21 


Total 

i 23 

59.81 




Significant at 0.95 level. 

★ 

Adjusted multiple correlation coefficient: 0.749. 


2-16 


















TABLE 2*4.— STATISTICS FOR MODEL 4 
(a) Regression coefficient 


Variable 

Coefficients* 

Standard error 

t-value 

1 

1.05758 

0.04131 

25.60* 

-n 



Constant 

3.82065 

5.30875 


X 1 

.01236 

.09462 

.13 

x 2 

-.16318 

.21757 

-.75 

x 3 

X 4 

.23395 

.48307 

.12736 

.19035 

1.84 

2.54 + 

4 

X 5 

x 6 

-.20947 

.59515 

.06269 

.07180 

-3.34 

8.2<i + 

tn (constant) 

.43214 

. 28467 



(b) Analysis of variance 


Source 

Degrees of freedom 

Sum of squares 

Mean of squares 

F-value 

Regression 

1 

1 1 26673 88 

1126673.88 

655.15 

Residual 

122 

209803.89 

1719.70 


Total 

123 

1336477.77 




Significant at 0.95 level. 

★ 

Adjusted correlation coefficient: 0.917. 
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potential of CMI and the mean pupation time were significant. The growth 
potential of the short-term mean air temperature was weakly significant. The 
growth potential of the long-term air temperature was not significant, prob- 
ably due to the high correlation (0.825) with the growth potential of mean 
pupation time. 

Model 2 adds the parent generation Infestation as an Independent variable to 
the set of weather parameters In Model 1. Model 2 also predicts Infestation. 
The adjusted multiple correlation coefficient was 0.708; the growth potentials 
of short- and long-term mean air temperatures were not significant In this 
model . 

Model 3 adds a sixth Independent variable In the form of the observed growth 
per generation 1 week prior to the prediction time. This model predicts 
growth per generation. The adjusted multiple correlation coefficient was 
0.749. As In Model 2, the growth potentials of short- and long-term mean 
air temperatures were not significant. 

The differences between Model 1 and Model 3 are important and should be 
clearly understood. Model 1 predicts the potential growth per generation 
as a function of weather. Thus, a high potential for growth may exist any- 
where if weather conditions are favorable. Model 3 predicts actual growth 
per generation of an existing population. If no screwworms are present in an 
area at the time of the parent generation, Model 3 will predict no growth. 

Model 4 fits the product of the growth predicted by Model 3 and the observed 
infestation during the parent generation with a simple linear regression to 
predict Infestation. The adjusted multiple correlation coefficient was 0.917. 
The growth potentials of the short- and long-term mean air temperatures 
remained nonsignificant due to the high intercorrelation of the three 
temperature-based terms. 

This model has an interesting feature in that, due to the constant term, it 
will introduce a small (three cases) infestation into an area regardless of 
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the weather conditions. However, tne population will never grow unless the 
condition becomes favorable. 

The accuracy of predicted infestations Is one measure of the relative effec- 
tiveness of these models. Table 2-5 gives the relative performance of each 
model as measured by the standard error of estimate. 

TABLE 2-5.- STANDARD ERROR OF ESTIMATE FOR EACH MODEL 


Model 

Cases 

Mean, % 

1 

79.1 

0.949 

2 

70.3 

.844 

3 

41.6 

.413 

4 

41 .4 

.411 


Models 2 and 4 are case models and can be evaluated directly. The predicted 
cases from the growth Models 1 and 3 were calculated by multiplying the 
predicted growth by the number of cases observed during the parent generation. 
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3. CONCLUSIONS 


Although complete validation of models of this type Is difficult and requires 
extensive field data. It Is clear that relatively simple weather-based models 
can successfully be used. A significant portion of the variation In screw- 
worm population growth and development has been traced to weather-related 
parameters. Clearly, valuable Insights Into the complex ecological situation 
governing screwworm population dynamics are available through simple weather 
models. 
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